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Ohmic contacts to a two-dimensional electron gas (2DEG) in GaAs/Ala;Gai_a;As heterostructures 
are often realized by annealing of AuGe/Ni/Au that is deposited on its surface. We studied how 
the quality of this type of ohmic contact depends on the annealing time and temperature, and 
how optimal parameters depend on the depth of the 2DEG below the surface. Combined with 
transmission electron microscopy and energy-dispersive X-ray spectrometry studies of the annealed 
contacts, our results allow for identifying the annealing mechanism and proposing a model that can 
predict optimal annealing parameters for a certain heterostructure. 



I. INTRODUCTION 

Epitaxially grown GaAs/Ala;Gai_2;As heterostructures 
that contain a two-dimensional electron gas (2DEG) 
are widely used for ele ctron transport studies in low- 
dimensional system|3^. Establishing electrical contacts 
to the 2DEG is a crucial step in device fabrication with 
these heterostructures. A commonly used recipe for mak- 
ing ohmic contacts is annealing of a AuGe/Ni/Au alloy 
that has been deposited on the heterostructure surface'^. 
High-quality heterostructures are often only available in 
a limited quantity, and it is desirable to minimize the 
heating that is needed for annealing the contacts to avoid 
damaging the heterostructure. A model that predicts 
optimal annealing times and temperatures for a het- 
erostructure with the 2DEG at a certain depth is there- 
fore very valuable. 

We present here a study of the annealing mechanism 
for this type of ohmic contact, and a model that can 
predict optimal annealing parameters for a certain het- 
erostructure. We used electron transport experiments to 
study how the quality of AuGe/Ni/Au based ohmic con- 
tacts depends on annealing time and temperature, and 
how the optimal parameters change with the depth of 
the 2DEG below the surface. These results confirm that 
the annealing mechanism cannot be described by a single 
simple diffusion process. Cross-sectional studies of an- 
nealed contacts with Transmission Electron Microscope 
(TEM) and Energy Dispersive X-ray (EDX) techniques 
were used for identifying a more complex annealing mech- 
anism, that is in agreement with the results from our 
electron transport studies. 

The AuGe/Ni/Au contact was first introduced by 
Braslau et al!^ to contact n-GaAs, and several stud- 
ies aimed at understandi ng the contact mechanism 
for this type of contaci! ^ l t' l 7m^m" | i^ l " | i4 | i5 | i l|i7|il^ 

Later studies focussed on the formation of an 
ohmic contact to a 2DEG in a GaAs/Al^jGai^ajAs 
heterostructure| i" l ^" l ^i | ^^l^^l^4 | 25 ]2ll27|^ but do not report 

how the optimal annealing parameters depend on the 



depth of the 2DEG below the surface. A number of these 
studies suggest that a contact is formed because a pat- 
tern of Au/Ni/Ge spikes that originate from the metal- 
lization penetrate the h eterostructure, just beyond the 
depth of the 2DEGP^l2ll, Earlier work had already sug- 
gested that in good contacts Ni-rich phases may form at 
the depth of the 2DEG, in contact with the GaAs below 
the Al^Gai„a;Aji2!, observe, instead, a mechanism 
where metal-rich phases only penetrate the heterostruc- 
ture over a distance that is shorter than the depth of the 
2DEG. The mechanism that results in a good contact is 
then similar to a process that has been describecP for 
contacts to n-GaAs: during annealing, the AuGe/Ni/Au 
on the surface segregates in Ni-rich and Au-rich domains, 
where the Ni domains contain most of the Ge. These do- 
mains penetrate the heterostructure and grow towards 
the 2DEG in large grains rather than narrow spikes. For 
optimal electrical contact conditions the Au and Ni-rich 
grains do not reach the 2DEG. The contact resistance 
decreases and the contact becomes ohmic because Ge dif- 
fuses deeper, forming a highly doped Alj;Gai_a;As region 
between the 2DEG layer and metal-rich phases at the 
surface. We find that even for very long annealing times, 
when the contact resistance has significantly increased 
compared to the optimal contact, the Au and Ni-rich 
phases still do not penetrate the 2DEG. 

The outline of this article is as follows: we first de- 
scribe our wafer materials and device fabrication. Next, 
we present electrical measurements and use these to iden- 
tify annealed contacts with optimal ohmic properties. In 
Section |IV| we present the results of our TEM and EDX 
studies of annealed contacts. Section |V] then summarizes 
the annealing mechanism that we identified, and this is 
used in Section |VI| to propose a model that can predict 
optimal annealing parameters. Finally, in Section [VII[ we 
present a study of how the contact resistance depends on 
the shape of the contact (varying area or circumference), 
which gives further insight in the annealing mechanism 
and the electrical contact properties. 
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II. DEVICE FABRICATION 

We studied annealed AuGe/Ni/Au contacts to three 
GaAs/Al2:Gai_a;As heterostructures, grown on (001)- 
oriented i-GaAs substrates, with the 2DEG at a hetero- 
junction at 70 nm (wafer A), 114 nni (wafer B), and 
180 nm (wafer C) below the surface of the wafer. These 
wafers have similar values for the 2DEG electron density 
Ug and mobility /j, (around 2 • 10^^ m^^ and 100 m^/Vs, 
respectively, results for 4.2 K and samples kept in the 
dark during cool down). For all three wafers the layer 
sequence (from the surface down) is very similar besides 
the depth of the 2DEG. The top layer is a ^ 5 nm n-GaAs 
capping layer, then an Alj^Gai-^As doping layer (Si at 
~ 1 • 10^^ cm-3) with X « 0.32, of thickness 30 nm (A), 
72 nm (B) or 140 nm (C). After this follows an undoped 
Al:rGai_^As buffer layer (~35 nm thick). The 2DEG is 
located at the interface with the next layer, which is a 
several fim thick undoped GaAs layer. 

We studied 200 x 200 /^m^ contacts that were de- 
fined by optical lithography on a 1 mm wide and 2 mm 
long etched mesa with a typical Hall-bar geometry. An 
electron-beam evaporator was used for deposition of sub- 
sequently 150 nm AuGe of eutectic composition (12 wt% 
Ge), 30 nm of Ni and 20 nm of Au. Subsequent anneal- 
ing took place in a pre-heated quartz furnace tube with 
a clean N2 flow of about 1 cm/s over the sample to pre- 
vent oxidation. We also found that using a much weaker 
N2 flow could result in surface contamination that was 
electrically conducting. We have used three different an- 
nealing temperatures, 400 °C, 450 °C and 500 °C. The 
samples were placed on a quartz boat and then moved 
into the center of the oven for various annealing times. 
Figure [iji shows the temperature of the surface of the 
quartz boat as a function of time for an oven tempera- 
ture of 450 °C. We assume that the sample temperature 
closely follows the temperature of the quartz boat, since 
we assured a good thermal contact over the full surface 
of the sample. 



III. ELECTRICAL MEASUREMENTS 

We measured the current-voltage (IV) characteristics 
of all contacts to determine optimal annealing parame- 
ters. We found that a suitable and sufficient definition 
for an optimal ohmic contact is a contact with the lowest 
zero-bias resistance at 4.2 K. The typical resistance for 
such a contact is ^--^20 il, but we have observed resistances 
as low as 5 ri. These values for contact resistance are 
close to the lowest values that have been reported—. All 
contacts defined as optimal in this manner showed highly 
linear IVs up to at least 1 mV (over and under annealed 
contacts did show non-linear IVs due to effects like Schot- 
tky or tunnel barriers in the contacts). Further, all these 
optimal contacts showed a strong monotonous reduction 
of the contact resistance upon lowering the sample tem- 
perature from 300 K to 4.2 K. Highly over and under 



annealed contacts showed an increase of the contact re- 
sistance upon cooling to 4.2 K. 
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FIG. 1: (a) Temperature of the quartz boat as a function of 
time for an oven temperature of 450 °C. Horizontal dashed 
line indicates the AuGe melting temperature T = 363 °C. The 
vertical dashed line indicates our definition of the annealing 
time tA, the time at which the boat is taken out of the oven, 
(b) Average contact resistance (i?) as a function of anneal- 
ing time Ia for contacts on wafer C, anneaied at 450 °C. A 
parabolic fit is made to estimate the annealing time where the 
resistance has a minimum, (c) Overview of optimal annealing 
times tA,Opt as a function of depth d of the 2DEG beneath 
the wafer surface for T = 400, 450 and 500 °C. The three gray 
lines (bottom to top for 400, 450 and 500 °C) are results of 
fitting a simple diffusion model to the experimental data (see 
text), which does not yield a good fit. 

We used a current-biased 4-terminal configuration to 
measure the voltage drop across a single contact, with 
two terminals attached to the metal bond wire that is 
pressed onto the ohmic contact, and two terminals at- 
tached to the 2DEG via other ohmic contacts (this al- 
lowed us to use a standard sample design in our fabri- 
cation facility). We accounted for a small voltage drop 
in the 2DEG area between the contact and the voltage 
probe via the 2DEG . We are aware that the Transmission 
Line Method (TLWpim 

is a better method for determin- 
ing the exact value of a contact resistance, but this is not 
needed for our approach. We compare resistances of var- 
ious annealed contacts that were fabricated under iden- 
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tical conditions besides the variation in anneaHng time 
and temperature. Within such a set, we determine which 
contacts have the lowest contact resistance. When repro- 
ducing our results with contacts that were fabricated in 
a different batch (using the same electron-beam evapo- 
rator, but after replenishing the AuGe target), we find 
that the values of the lowest contact resistance can be 
different up to a factor 2 around the typical result. We 
attribute these batch-to-batch fluctuations to variations 
in the exact composition of the AuGe/Ni/Au layer that 
we deposit. The optimal annealing times, however, show 
batch-to-batch fluctuations of only 10%. Thus, our ap- 
proach for determining optimal annealing conditions does 
not depend on the exact value of the measured contact 
resistance. 

Figure [T]d shows a typical result, from which we deter- 
mine the optimal annealing time for contacts to wafer C 
for the case of annealing with the oven at 450 °C. Contact 
resistance data that is denoted as (R) is the average resis- 
tance measured on a set of 8 identical contacts, and the 
error bar represents the standard deviation. The results 
in Fig. show a clear minimum in contact resistance 
for annealing times near 5 minutes. We fit a parabola 
(phenomenological ansatz) to the log(i?) values of these 
data points, and define the optimal annealing time as the 
time coordinate of the minimum of the parabola. In this 
manner, the optimal annealing times tA,opt are obtained 
for contacts on wafers A, B and C, annealed at each of 
the temperatures. 

Figure [TJ; presents these optimal annealing times. As 
expected, the optimal annealing time increases as the 
temperature is decreased, and increases as the depth d 
of the 2DEG increases. While it is known that several 
simultaneous diffusion processes play a role in contact 
formation^, we will, for the sake of argument, show that a 
simple diffusion model has little value for predicting how 
optimal annealing times depend on the depth d and the 
annealing temperature. For this simple diffusion model, 
we assume that a certain dopant (with fixed concentra- 
tion Co at the surface) diffuses into the heterostructure. 
The relevant solution to Fick's second law is then 



C = Co erfc- 



(1) 



Here C is the doping concentration at time t and depth z 
into the heterostructure, and D is the diffusion constant 
(erfc is complementary error function). Since the tem- 
perature of our sample is not constant (see Fig. [iji) we 
will use the measured temperature profile T(t) to inte- 
grate the diffusion constant over time, and use in Eq. [l] 
J D{t)dt instead of Dt, where 



D{t) = Do exp(- 



Ea 



ksTit) 



(2) 



and where Ea is an activation energy. We assume that 
an optimal contact then always occurs for a certain op- 
timal value for C/Cq at the depth of the 2DEG {z = d). 
We define the annealing time as the time from start to 



the moment when the boat is taken out of the oven, but 
integrate over the entire time span that the sample is 
at elevated temperatures, (as shown in Fig. [iji, fully in- 
cluding the cooling down). This gives a model with the 
activation energy Ea, diffusion constant Do and concen- 
tration C /Co as fitting parameters. 

The gray lines in Fig. [l]; show the best fitting result 
that reasonably covers all 9 data points in a single fit. 
Besides the fact that the shape of the traces only poorly 
matches the trend in the data, the parameter values give 
unreasonable results. The temperature dependence alone 
governs the fitting result for Ea, giving here 0.15 eV. This 
is on the low side for typic al values for diffusion in GaAs 
materials (^1 eV)EiE2l33l^ por fixed Ea, various combi- 
nations of C/Cq and Do give identical results. When 
assuming a typical value D o ^ 3 ■ 1 0~^ m^/s (for diffu- 
sion of Ge, Ni or Au in GaAspHMM]) ^ ^j^jg yjgjjg (jjfj^ 

very close to 1, i.e. completely saturated diffusion. This 
is in contradiction with the clear dependence on depth 
that we observe (and this remains the case when allow- 
ing for Ea up to ~ 1 eV, but then the fit does not cover 
all 9 data points at all). Thus, we find that predicting 
optimal annealing times with simple diffusion (according 
to tA,Opt oc d^ at constant temperature) does not work 
and that a more complex model needs to be considered. 



IV. TEM AND EDX RESULTS 

We have studied the contact formation using cross- 
sectional TEM imaging of contacts at several stages dur- 
ing the annealing process. The samples were prepared 
for TEM imaging by using a Focussed Ion Beam (FIB) 
to slice out a micrometer thin piece of the measured con- 
tact. By further thinning using the FIB the thickness 
was reduced to 100 nm. 

Figure |2^ shows an overview of an optimally annealed 
contact on wafer C which was annealed for 5 minutes at 
450 °C. The composition of the various phases has been 
determined by Energy Dispersive X-ray (EDX) analysis 
and is illustrated in Fig. ^p. From bottom to top we 
recognize the GaAs substrate, and an AlAs/GaAs su- 
perlattice to smoothen the surface of the substrate. On 
top of that we find another layer of epitaxially grown 
GaAs and a layer of Ala;Gai_a;As. The 2DEG is at the 
interface of these two layers. The GaAs capping layer 
that was originally on top of the Ala;Gai_a;As layer is no 
longer visible. Instead, we see large Au-rich and Ni-rich 
grains that have penetrated below the original wafer sur- 
face. Both of these phases contain out-diffused Ga and 
As, with Ga mainly in the Au-rich grains and As mainly 
in the Ni-rich grains. Further, the Ni-rich phase absorbed 
most of the Ge. We find that the Au grains do not con- 
tain any Ge, consistent with the findings of Kuan et alW 
in work with n-GaAs. 

The wide and curved dark lines going over all the het- 
erostructure layers (most clearly visible in the GaAs lay- 
ers) are due to strain induced by the FIB sample prepa- 
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FIG. 2: (a) Cross-section TEM image of a contact on wafer C, 
annealed for the optimal annealing time at 450 °C. (b) A 
sketch of the TEM image in (a) to specify the various layers 
and phases, (c) Larger area TEM image of the same contact 
as in (a) showing large Au-rich (black) and Ni-rich grains 
(dark gray) contacting the Al3;Gai_a;As. (d) Similar image 
for a highly over annealed contact. The Au and Ni grains still 
do not penetrate the 2DEG, but Au has diffused underneath 
the Ni grains, which results in an increased contact resistance. 



ration and are not related to the diffusion process. 

We find that the Au-rich and Ni-rich grains do not 
have to penetrate the 2DEG in order to establish a good 
electrical contact. We can rule out that we do not see 
grains reaching the 2DEG due to the small thickness of 
the sample slice, since we observed no substantial vari- 
ation in the penetration depth of a large number of Au 
and Ni grains going along the sample slice. We examined 
two slices from two different samples, both with a length 
of 100 /im, after electrical measurements confirmed that 
these contacts were indeed optimally annealed. 

The TEM image in Fig. [2]: shows a larger region of an 
optimally annealed contact. Large Au and Ni grains that 
have penetrated the Al2;Gai_a;As layer can be identified. 
Figure [2}i shows an over annealed contact on wafer C, 
that was annealed for 7 minutes at 450 °C. Remarkably, 
the Au and Ni grains did not penetrate much further into 
the Al2;Gai_a;As, and do still not reach the 2DEG. The 
most significant change with respect to Fig. [2]: is that the 
Au-rich phase is diffusing underneath the Ni-rich grains, 
reducing the total Ni-grain-Ala;Gai_3;As interface area. 
This was also observed by Kuan et al^ (and confirmed in 
detailed studies by Lumpkin et alJ^ in work on n-GaAs, 
and the results of these authors indicate that this process 
is mainly responsible for the increase in contact resistance 
when a sample is being over annealed. 



Kuan et aZ.l^l report that the contact resistance is sen- 
sitive to the ratio of the total contact area between Au- 
rich regions and Al^Gai_a:As, and that of Ni-rich regions. 
The Au-Alj;Gai_2:As interface is considered a region of 
poor conductance because the Au-rich grains (in contrast 
to Ni-rich grains) do not contain any Ge, such that it 
cannot act as a source for diffusion of Ge into the het- 
erostructure. However, it is to our knowledge not yet 
understood why the diffusion of Au underneath the Ni 
grains at later stages of annealing (when a large amount 
of Ge already diffused out of Ni) results in a strong in- 
crease of the contact resistance. 



V. SUMMARY OF ANNEALING MECHANISM 

In this Section we use the results from the previous two 
Sections, together with established results from the liter- 
ature, for giving a qualitative description of the formation 
of an ohmic contact to a 2DEG in a GaAs/Ala;Gai_2;As 
heterostructure. It is remarkably similar to the anneal- 
ing mechanism as described by Kuan et al^ for contacts 
to n-GaAs. In the initial stages of the annealing pro- 
cess (already during the heating of the sample) Au and 
Ge segregate, and most Ge forms a new phase with the 
Ni. At the same time, these Ge-rich Ni grains move to 
the wafer surface due to a wetting effecifiS, which results 
in the situation that the wafer surface is covered with 
neighboring Au and NixG^^ grains. There is evidence 
that for thin metallization layers ('^lOO nm) this process 
already occurs well below the bulk melting temperature 
(363 °C) of the eutectic AuGe phastP. 

Next, at higher temperatures, both the Au-rich and 
Ni-rich grains penetrate into the heterostructure by solid 
phase inter diffusion, compensated by a back flow of As 
and Ga. Our EDX results confirm that Ga mainly flows 
into Au, and As mainly into Ni-rich grains. This concerns 
the forma tion of new ma terial phases. In several ear- 
lier studiej^ ^^l"'^^!^^ ! -'^'' ! ^^!-'^^! these phases have been identi- 
fied as AuGa alloys and phases close to Ni2GeAs. These 
phases penetrate only tens of nm below the original wafer 
surface for typical annealing conditiond^El^ 

At the same time, there is diffusion of atomic Ge, Ni 
and Au (at similar concentrat ions) into t he heterostruc- 
ture, which penetrates deepeilSESIiSMIIZl jn particular, 
Ge diffuses out of the Ni-rich grains into the Ala;Gai_a;As 
layers, and a good ohmic contact is formed when the 
AljjGai-jjAs layers are sufficiently doped with Ge all 
the way up to the 2DEG. While this is progressing, the 
Au-ri ch g rains start to expand underneath the Ni-rich 
graingSES, which have the lowest contact resistance with 
the doped Ala;Gai_a:As layer since they were the domi- 
nant supplier of Ge. The expansion of the AuGa grains 
is possibly due to the relatively low activation energy for 
out diffusion of Ga into AiP^l (while the Al-As binding 
energy is relatively high^^. This latter process increases 
the interface resistance between the metallization on the 
surface and the doped Al2;Gai_2:As layer. Thus, the for- 
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mation of an optimal contact is a competition between 
these two processes. 

The in-difFusion of Ge lowers the contact resistance for 
two reasons. 1) The full Alj.Ga1_2.As region between the 
surface and the 2DEG becomes a highly doped region 
with a reasonably low bulk resistivity. 2) The Ge dop- 
ing in this region makes the Schottky barrier between 
the doped semiconductor and the surface metallization 
very thin (the barrier height is probably not changing 
significantlj^^ , up to the point where its series contri- 
bution to the contact resistance is small. The total con- 
tact resistance is then dominated by doped Ala;Gai_a;As 
region, giving linear transport characteristics (a similar 
effect occurs for contacts to n-GaAs due to spreading re- 
sistance below the contact'^. 

As said, it is not yet well established which processes 
are responsible for the resistance increase upon over an- 
nealing. The fact that over annealing with 2DEG samples 
and n-GaAs samples^ occurs qualitatively in a very simi- 
lar manner (and also at similar annealing times and tem- 
peratures) is a first indication that it is due to a process 
near the interface with metal-rich phases on the surface, 
rather than a process at the depth of the 2DEG or the 
edge of a contact. Further, our results now show that 
the resistance increase for 2DEG samples is also corre- 
lated with the expanding AuGa grains below the Ni-rich 
grains. Various authors have suggested that the increas- 
ing contact resistance that is associated with over anneal- 
ing may be due to a large number of vacancies just below 
the metal-rich phases near the surface^^ (but others 
suggested it was due to excessive in-difi'usion of NP^ll). 
These mainly result from out-diffusion of Ga into the Au- 
rich grains (which indeed results in a very stable AuGa 
phase near the original wafer surfac^ ^^ * ^^ !). These va- 
cancies occur in particular when there is no (longer) Ge 
diffusion into these vacancies. One should note, however, 
that with n-GaAs an increasing contact resistance was 
also observed without an expansion of the AuGa grains 
below the Ni-rich grain^, but this does not rule out that 
an increasing number of vacancies is responsible for over 
annealing. 

Finally we remark that both the Ni-rich and Au-rich 
grains are probably important for rapid annealing at rel- 
atively low temperatures. The Ni-rich grains act as the 
supplier of Ge. The presence of Au grains may be im- 
portant since it rapidly results in a large number of Ga 
vacancies. This probably enhances the in-diffusion of Ge. 
It was for example also observed that the creation of such 
vacancies near the surface, enhances the diffusion of Si 
dopants from the doping layer (much deeper into the ma- 
terial) into neighboring layers^. 



VI. DIFFUSION MODEL 

We use the above description to construct a model that 
predicts the optimal annealing time for a given anneal- 
ing temperature and 2DEG depth d. The contact re- 



sistance is then the series resistance of the resistance of 
the Ge-doped Ala;Gai_a;As region (i?Ge) and the inter- 
face resistance between the surface metallization and this 
Ge-doped Ala;Gai_j;As layer (i?i/). For both we consider 
the average over the full contact area. We will first as- 
sume an anneal temperature T that is constant in time. 
We model the resistance of the Ge-doped AlxGai_j;As 
region using the result from work on n-GaAs that the 
contact resistance is inversely proportional to the doping 
concentration"^. Thus, we assume that 

where C{z)/Cq is the local Ge concentration at depth z 
as in Eq. [T] and where the integral runs from the depth 
of the Au and Ni grains to the depth of the 2DEG. The 
behavior of this equation is that Rqc first rapidly de- 
creases, and then curves off to saturate at a level that is 
proportional to d (dashed curve in Fig.|3^). 

To model Rif, we assume that the increase in resis- 
tance for over annealed contacts is related to the de- 
crease in Ni-grain-Al2;Gai_a;As interface area. Imag- 
ine, for simplicity, a single, square shaped Ni-rich grain 
with area A^i — Lj^^. We model the reduction of 
this area as a sideways diffusion process of Au, again 
with a time-dependence as simple diffusion analogues 
to Eq. [1] The length of a side is then reduced as 
LNi{t) ~ Lq ~ '^y/^DAut, where Lq is the initial grain 
size, and Dau the diffusion constant for this process, such 
that 

^ (Lo- 274^)2- 

For a very wide parameter range, this model gives that 
Rif increases more or less linearly in time (solid black 
curve in Fig. |3^). A resistance increase that is much 
stronger than linear only sets in when the total inter- 
face area approaches zero, when the contact is already 
strongly over annealed. The total contact resistance is 
the sum of i?(3e and Rif (gray solid curve Fig. ^), and 
the optimal annealing time is then defined as the time 
where this sum shows a minimum value. 

We can reduce the number of fitting parameters for 
this modeling to only two with the following approach. 
For i?Ge in Eq. [3] we assume parameters where Rce sat- 
urates at a value below, but on the order of the optimal 
contact resistance Ropt- We also assume that this satura- 
tion occurs in a time scale on the order of a few times the 
optimal annealing time. For Rif in Eq.ffl we assume that 
it has a value below Ropt for t = 0, ana that it increases 
more or less in a linear fashion to a value of order Ropt ■ 
This increase should take place in a time scale on the 
order of the optimal annealing time. Numerically inves- 
tigating this model then shows that it has for a very wide 
parameter range the behavior that the increase of opti- 
mal annealing time ^A.Opt with increasing 2DEG depth d 
is close to linear. We can express this using an effective 
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FIG. 3: (a) Model for the resistance of an ohmic contact as 
a function of annealing time at constant temperature. The 
resistance Roe of the Al^Gai_^As layers (dashed line) de- 
creases in time due to increased Ge doping. The interface 
resistance Rif between the surface metallization and the Ge- 
doped Ala;Gai_a;As layers (solid black line) increases in time 
due to a decreasing Ni-grain-AUGai-xAs interface area. The 
time where the sum of these two resistances (gray solid line) 
shows a minimum defines the optimum annealing time tA,Opt- 
(b) Effective velocity of optimal contact formation Vocf as a 
function of temperature (Eq. [6|, plotted for parameters that 
give the best fit in (c). (c) Optimal annealing times as the 
2DEG depth and annealing temperature is varied (same ex- 
perimental data as in Fig.|lJ;). The solid gray lines (left to 
right for 500, 450 and 400 G) represent fits using the model 
of Eqs. ^and[6] (see text for details). 



velocity for optimal contact formation Vocf, 

tA,Opt = d/Vocf- (5) 

Furthermore, numerical investigation of the temperature 
dependence shows that Vocf behaves according to 

Vocf{T) = vo exp(-— ^) (6) 

when the diffusion processes that underlie Eq.[3]and Eq.|4] 
are both thermally activated with a similar activation 
energy Ea- We can now fit this model to our experimental 
data only using Eq.[5]and Eq.[6] such that we only have vq 
and Ea as fitting parameters. In doing so, we take again 



into account that the temperature T{t) is not constant 
during annealing, and use again profiles as in Fig Ilk. 

The results of this fitting are presented in Fig. [3p, and 
Vocf as a function of temperature for these fitting param- 
eters {Ea — 0.6 eV and vq — 7.6 • 10~^ m/s) is plotted 
in Fig. |3]d. While it is a crude model, the fits are very 
reasonable, showing that the model is useful for predict- 
ing optimal annealing t imes. F urthermore, the value for 
Ea is a realistic numbeil^SIIIll] q^j. j^odel also predicts 
that the minimum value of the resistance that can be 
achieved for optimally annealed contacts increases with 
increasing 2DEG depth. We did not observe such a clear 
trend, probably because the resistance of optimal con- 
tacts is so low that one needs to include contributions 
from 2DEG square resistance around and underneath 
the contact when evaluating absolute values (further dis- 
cussed below). 



VII. CONTACT-SHAPE DEPENDENCE 

Our model for the annealing mechanism implies that 
optimal contacts have a rather uniform Ge concentration 
throughout the Al^jGai.j-As layers, and that this results 
in a value for of order 10 fi. This implies that the 
bulk resistivity in the doped Ge-doped Ala;Gai_2;As layer 
is around 4 fim. In turn, this implies that in-plane elec- 
tron transport under an optimal contact from the metal- 
lization on the surface to 2DEG on the side of the con- 
tact still mainly takes place in the original 2DEG layer. 
If the square resistance i?n for transport in the original 
2DEG layer below the contact does not strongly increase 
during annealing, and if it is smaller than the contact re- 
sistance, this also implies that the resistance of optimal 
contacts should be inversely proportional to the contact 
area. Thus, measuring whether the contact resistance 
depends on contact area or on the circumference of a 
contact can give further insight in the annealing mecha- 
nism and contact properties. 

We carried out such a study, by varying the shape 
of contacts. All results that we discussed up to here 
were obtained with square contacts with an area A of 
0.04 mm^ and a circumference Cl = 4L of 0.8 mm (on 
the side of a Hall bar). For the dependence on con- 
tact shape, we measured various sets where we varied 
the circumference Cl while keeping the area constant at 
0.04 mm^, and various sets where we varied the area while 
keeping the circumference constant at 0.8 mm. We varied 
the shape from smooth circular shape to square shapes 
with a zig-zag edge at the 50 micronscale, to avoid getting 
too much resistance contribution from square resistance 
of 2DEG right next to a contact (for these devices we 
used electron-beam lithography). The study only used 
wafer A. All contacts were fabricated and annealed in 
one single batch to ensure that it is meaningful to com- 
pare the values of contact resistance. 

For this study, we inject again current into the con- 
tact that is measured, and extract the current using 
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another contact. However, the dependence on contact 
shape can only give an unambiguous result if the resis- 
tance from each side of the studied contact to the place in 
the 2DEG where the current is extracted is sufficiently 
similar. This can be achieved by making the distance 
between the contacts larger than the size of the con- 
tacts. Thus, we now fabricated contacts in the mid- 
dle of 2 mm by 3 mm cleaved wafer pieces (two rows 
of four contacts, with center-to-center distance between 
rows 1 mm, and center-to-center distance between con- 
tacts within a row 0.6 mm). Using four different con- 
tacts for a 4-terminal measurement on the 2DEG (with 
the current biased from one row to the other) gives on 
these samples indeed low values around 8 fi, in reason- 
able agreement with the value of the 2DEG square resis- 
tance i?n of about 20 il. Contact resistance values were 
again determined in a current-biased 4-terminal configu- 
ration, with two terminals connected to the bond wire on 
the contact, the second current terminal on the opposite 
contact in the other row, and the second voltage terminal 
on a neighboring contact in the same row. 

On contacts that are not annealed, we can observe a 
tunnel current, as expected for Schottky barriers. Here, 
the effective resistance is inversely proportional to area. 
For optimally annealed contacts, we found that the con- 
tact resistance was independent on circumference, while 
only showing a weak dependence on area (weaker than 
inversely proportional to area), see Fig.|4] The fact that 
the dependence on shape does here not show a clear de- 
pendence as (i?) oc 1/A agrees with the fact that the 
(R) values are comparable to the square resistance of 
the 2DEG, such that the latter gives a significant con- 
tribution to the total contact resistance. Fully under- 
standing the contact resistance then requires incorporat- 
ing all square resistance contributions from underneath 
and around the 2DEG. Since we found it impossible to 
estimate these effects with a small error bar, we tried to 
demonstrate a clear dependence on area by measuring 
sHghtly under annealed contacts instead. 
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FIG. 4: Contact resistance (R) as a function of (a) contact 
area A for constant circumference 41/ and (b) contact circum- 
ference C for constant area A. The error bars here represent 
the standard deviation from measuring i? on 8 identical con- 
tacts. The dashed line in (a) is a fit using (R) oc 1/A. 

On two sets of under annealed contacts on wafer A, 
where we used shorter anneal times than tA,opt (average 



contact resistance of 30 and 500 fi), we found (within 
error bar) no dependence on area or circumference. We 
can only explain this result if we assume that the 2DEG 
square resistance underneath the contact is significantly 
increased (to values comparable to the total observed 
contact resistance) for under annealed contacts. This 
probably results from the in-diffusing Ge (and atomic Au 
and Ni^^ ■^^), which already introduces strain and scatter 
centers in the 2DEG layer before optimal contact condi- 
tions are reached. For optimal annealed contacts (here 
with total resistance of typically 7 fl, independent of cir- 
cumference), the square resistance underneath a contact 
must have returned to a low value of order 10 f2. Appar- 
ently, the resistance increase due to strain and scatter 
centers is compensated by increased Ge doping near the 
2DEG layer. 

The summary of this study is that the resistance of 
annealed contacts never shows in a clear dependence on 
circumference, and only a weak dependence on area for 
optimal contacts. We can, nevertheless, draw the follow- 
ing conclusions. For an optimal ohmic contact, it is not 
the case that electron transport between the surface met- 
allization and the surrounding 2DEG mainly occurs at 
the edge of a contact. Instead, the full contact area plays 
a role, and in-plane electron transport under an optimal 
contact mainly takes place in the original plane of the 
2DEG. In addition, we find it impossible to evaluate the 
absolute values of the contact resistance of our devices 
with an accuracy within a factor 2, since the resistance 
of an optimal contact has a contribution from the square 
resistance underneath the contact, and its value is influ- 
enced by the annealing process. Further, future studies in 
this direction should consider that pressing the bond wire 
(with a footprint of about 100 x 100 /im^) onto the surface 
metallization may locally disturb the contact properties, 
which can disturb a clear dependence on contact shape. 

We could therefore not study the property of our model 
that the optimal contact resistance value should be pro- 
portional to d. Instead, we should evaluate whether the 
enhanced square resistance underneath a contact needs 
to be incorporated in our model. We find that this is 
not needed for the following reasons: for over anneal- 
ing it does not play a role, since we observe the same 
over- annealing mechanism as observed on bulk n-GaAs. 
Optimally annealed contacts occur when the square re- 
sistance underneath the contacts has again low values of 
order 10 f2. Here we observe a weak area dependence 
and no dependence on circumference, such that we can 
rule out that the effect dominates the contact resistance. 
Thus, the only effect is that it temporarily enhances the 
total contact resistance by about a factor 2 while the an- 
nealing progresses towards optimal contact conditions. 
Note that this does not change the fact that lowering the 
contact resistance in this phase still fully depends on fur- 
ther Ge diffusion towards the 2DEG layer. Therefore, it 
only slightly modifies how Rq^ in Eq.[3]decreases towards 
low values. 
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VIII. CONCLUSIONS 

Summarizing, we have measured the zero-bias resis- 
tance of annealed AuGe/Ni/Au ohmic contacts to a 
2DEG as a function of annealing time and temperature. 
We have thus obtained optimal annealing parameters for 
three different heterostructures where the 2DEG lies at 
a different depth below the wafer surface. TEM images 
of several annealed contacts provided further insight into 
the annealing mechanism and the formation of a good 
ohmic contact. 

Combining this information we have developed a 
model that can predict the optimal annealing param- 
eters for contacting a 2DEG at a certain depth in a 
GaAs/Al^^Gai-ajAs heterostructure. The model assumes 
two competing processes: 1) Diffusion of Ge into the het- 
erostructure lowers the contact resistance, and results in 
linear transport characteristics. 2) At longer annealing 
times, Au-rich phases diffuse in between the heterostruc- 
ture and Ni-rich phases at the wafer surface. The associ- 
ated increase in contact resistance is probably due to sub- 
sequent diffusion of Ga into this Au-rich phase, since this 
increases the number of Ga vacancies in the heterostruc- 
ture near the metallization on the surface. The competi- 
tion between these two processes results in a mechanism 
where the optimal annealing time (for a process at con- 
stant annealing temperature) is proportional to the depth 
of the 2DEG below the surface, and the speed of this pro- 
cess has thermally activated behavior. This model should 
have predictive power for many heterostructures, as long 
as the temperature of the samples as a function of time 



during the annealing process is known. Our study of 
how the contact resistance depends on the shape of the 
contact confirmed that the full contact area plays a role 
in electron transport between the metallization on the 
surface and the 2DEG. 

Our model may become invalid for systems with a very 
deep 2DEG, since Rif (Eq. [2]) is expected to increase 
more rapidly at long annealing times, possibly result- 
ing in non-ohmic behavior. Our results suggests that 
for solving this problem the focus should be at main- 
taining sufficient contact area between Ge-rich Ni grains 
and the Ge-doped Ala;Gai_2:As layer at long anneahng 
times. This can possibly be engineered by changing 
the layer thickness, order an d composition of the ini- 
tial AuGe/Ni/Au metallizatioiJiSlIlIillMiMmil^ or by in- 
cluding a Pt, Nb or Ag layer below the top Au layer that 
suppresses t he intermixi ng of this Au with layers at the 
wafer surfac^^E^E^'^^'21 (uniform Ni/Ge/As layers been 
reportecPil) . Alternatively, one can reduce the depth of 
the 2DEG by etching before deposition of AuGe/Ni/Au 
(up to the point where this results in depletion of the 
2DEG). 
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